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Abstract 

A toy model is proposed for four dimensional non-abelian gauge theories coupled to 
a large number of fermionic degrees of freedom. As the number of flavors is varied the 
gauge theory may be confining, walking or conformal. The toy model mimicking this 
feature is the two dimensional 0(3) sigma model with a i?-term. For all t? the model is 
asymptotically free. For small $ the model is confining in the infra red, for i? = n the 
model has a non-trivial infra red fixed point and consequently for ■& slightly below n the 
coupling walks. The first step in investigating the notoriously difficult systematic effects of 
the gauge theory in the toy model is to establish non-perturbatively that the ■& parameter 
is actually a relevant coupling. This is done by showing that there exist quantities that are 
entirely given by the total topological charge and are well defined in the continuum limit 
and are non-zero, despite the fact that the topological susceptibility is divergent. More 
precisely it is established that the differences of connected correlation functions of the 
topological charge (the cumulants) are finite and non-zero and consequently there is only 
a single divergent parameter in Z($) but otherwise it is finite. This divergent constant can 
be removed by an appropriate counter term rendering the theory completely finite even at 
■d > 0. 
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1 Introduction 



Lattice simulations of technicolor inspired models are plagued by known systematic uncer- 
tainties [TH5]. Although the models under consideration are QCD-like in that they are four 
dimensional non-abelian gauge theories coupled to dynamical fermions the systematic effects 
of the interesting models (those that are either conformal or walking) are much more difficult 
to control than in actual QCD. As a result currently there are disagreements between various 
approaches, discretizations, etc, and universality is not immediately evident [BT-llO]. Clearly the 
general expectation is that once all systematic effects are controlled and taken into account the 
results from different approaches and regularizations will agree as they should. 

In this paper a toy model is proposed which mimics many of the features of non-abelian 
gauge theories in the hope that systematic effects can be fully explored. Hopefully these will 
help controlling the corresponding effects in the much more complicated gauge theories. The 
proposed model is the two dimensional 0(3) non- linear sigma model with a d term. At *d = 
the model served as a toy model of QCD for a long time since it is asymptotically free, features 
instantons, confinement and dimensional transmutation [11) . It is exactly solvable |12j even at 
finite volume [13H15] . Since the topological term is invisible in perturbation theory the model 
is asymptotically free for arbitrary i). The dynamics in the infra red is however expected to be 
very sensitive to 

At ~ 7r the model is conjectured [161117] to have a non-trivial infra red fixed point governed 
by the SU(2) WZNW model at level k = 1 and, if the conjecture holds, is also exactly solvable. 
Some numerical evidence in support of the conjecture has been presented in [18] and a recent 
very detailed study confirming it in [19] . The infra red fixed point implies a zero of the j3- 
function. This situation is analogous to gauge theories in the conformal window. 

For < $ < ir exact solvability is lost but based on continuity one expects that for $ not 
much below tt the /^-function develops a near zero and the renormalized coupling will walk. This 
arrangement is analogous to gauge theories just below the conformal window. Hence dialing -d 
corresponds to dialing the number of flavors Nt in the gauge theory 

In all three scenarios (confining, walking, conformal) one may also introduce an external 
magnetic field to mimic the effect of a finite quark mass. 

Before exploring the analogies further and investigating the origins of the severe systematic 
effects the first task is to establish non-perturbatively that the $-term is actually a relevant 
operator and also what the singularity structure of the theory is for & > 0. This is not immedi- 
ately obvious largely because of the unusual scaling properties of the topological susceptibility 
and a class of similar observables. 

It is well known that small size instantons render the topological susceptibility \ = (Q 2 ) /V ill 
defined in the semi-classical approximation [20] . Going beyond the semi-classical approximation 
fully non-perturbative lattice studies have shown that regardless how one improves the details 
of the lattice implementation a logarithmically divergent susceptibility is obtained at finite 
physical volume in the continuum limit. Moreover, all even moments of the total topological 
charge distribution (Q 2m ) /V have the same property. 

However, the model at $ = is exactly solvable and both the exact solution and the 
continuum limit of lattice simulations agree that correlators of the topological charge density, 
e.g. (q(x)q(0)) are finite. The above two observations, namely that certain statistical properties 
of the total charge distribution P(Q) are ill defined while at the same time correlators of q(x) 
are finite, might make one wonder whether the total charge operator Q is an irrelevant operator 
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while q(x) is not. If so, the only consistent continuum value of (Q 2m ) would be zero and the 
apparent divergences in the lattice calculations would be regarded as artifacts. This scenario 
would imply that the theory defined on the lattice at non-zero ■d leads to an identical continuum 
theory as the one defined at $ = 0. Equivalently, the total charge operator inserted into any 
correlation function would be zero in the continuum theory (Q . . .) = 0, while correlation 
functions of the type (q(x) . . .) are finite. This scenario would of course invalidate Haldane's 
conjecture about the equivalence of the i? = ir theory with a non-trivial interacting conformal 
field theory. 

In this work it is shown that there exist quantities built out of the total topological charge 
operator Q which have well defined continuum limits and are non-zero. These observables 
are differences of connected correlation functions of the topological charge, in other words the 
cumulants. Each term is logarithmically divergent but the divergence cancels in the difference 
and moreover they scale correctly in the continuum limit to non-zero values. Showing correct 
scaling towards the continuum limit in itself would not be sufficient to prove that the i9-term 
is a relevant operator because the continuum limit value could be zero. Since all cumulant 
differences are finite there is only a single UV-divergent parameter in the partition function 
Z(€) but otherwise it is finite. 

While preparing this manuscript the preprint [19] appeared also with the conclusion that $ 
is a relevant coupling. The method was different though, in |19j it was shown to high precision 
that a well defined observable is different in the continuum limit for three different values of 
$ implying that $ can not be irrelevant. In the current work all simulations are carried out 
at = and the same conclusion is reached by showing that certain combinations of the 
topological charge operator are non-zero in the continuum. 

2 0(3) sigma model with a i9-term 

The model in Euclidean continuum notation is defined by the action 



for the unit 3- vectors s, s\ + s\ + s§ = 1, where go is the bare coupling. Only a torus geometry 
will be considered corresponding to a box of finite linear size L which will be regularized by a 
symmetric lattice. 

The corresponding partition function, free energy per unit volume and topological charge 
distribution of the model at non-zero "d and volume V is given by 



with the normalization Z(0) = P{Q) = 1- Since physics is periodic with period 2tt in $ 
and $ — > is a symmetry the free energy per unit volume can be Fourier expanded 




(1) 




(2) 



Q 



OO 




(3) 
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It has been pointed out in [3T] that in the semi-classical or dilute gas approximation all /„ 
coefficients vanish except for f\ which is UV divergent due to instantons of size a <C p <C £ where 
a is the lattice cut-off and £ is the physical correlation length. The remaining coefficients come 
from interactions between instantons. Semi-classical arguments also suggest that for instantons 
causing the UV divergence in /i the ratio between their size and their average separation goes 
to zero in the continuum limit. This would imply that the interactions responsible for the f n >i 
coefficients are small in the continuum limit hence will not cause them to diverge. 

To summarize, the semi-classical approximation accounts for a UV divergent f± and finite 
fn>i coefficients. A suitable way of addressing whether this statement is true beyond the semi- 
classical approximation is to consider observables that can be expressed by the f n >i coefficients 
only and calculating them fully non-perturbatively. The simplest choice is to take the connected 
correlation functions of the topological charge, 



X2m = {-If 

and consider their differences, 



d 2m f 

(4) 



(it) 



2m 



= 



Ax2m = X2m - X2m+2 = ^T? m {r? - 1) (5) 



n=2 

from which fx drops out. The first few such correlation functions are 

(Q 2 ) 



X 2 



V 

(Q 4 )-3(Q 2 ) 2 



X4 = ^ ' v ^ ' (6) 
(Q 6 )-15(Q 4 )(Q 2 )+30(Q 2 ) 3 

X6 = y ■ (7) 

All of these are expected to diverge in the continuum limit but their differences are expected 
to be finite. Some numerical evidence has been presented in [^TJ in favor of correct scaling 
behavior for A%2 but whether the continuum value is zero or non-zero has not been discussed. 

In the following it will be shown to high precision that the expectations from the semi- 
classical analysis indeed hold non-perturbatively and all moments (Q 2m ) and all cumulants 
X2m are logarithmically divergent but the differences A%2m are finite. This implies that there 
is a single ill-defined constant in F(d) namely /i but otherwise it is finite. The constant /i can 
be removed by an appropriate renormalization condition leading to a finite and universal free 
energy and partition function for arbitrary 



3 Numerical simulation 

It is convenient to take the continuum limit on a symmetric periodic lattice L 2 of fixed physical 
volume. Physical length and mass is defined by the second moment correlation length £2, 
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where 

M 2n = (f) E( 2sin r) ^ 



2^ 



is given in terms of the zero spatial momentum projection of the 2-point correlation function 
C(t,t') of the field s. Let us introduce m(L) = l/^iL). Note that in this notation m{L) is not 
the mass gap in finite volume but rather is simply defined as the inverse of £2 (which for L — > 00 
agrees with the mass gap but not for finite L). The physical volume is fixed to m(L)L — 4. A 
recently proposed [55] topological lattice action is used for the simulations, 

S = J2S(si, Sj ) (10) 

(i,3) 

where the sum is over all neighboring sites and 

Sis- s ■) = { ° if Sl ' Sj >C ° s6 (11) 
v *' 31 y 00 otherwise 

In other words the action is zero for two neighboring vectors if their relative angle is smaller than 
8 and infinite otherwise. The continuum limit is taken by tuning the bare coupling 8 towards 
zero. This action is topological because small perturbations of the field s do not change the 
action nevertheless it has been shown that it is in the right universality class |22j . 

If 8 < ir/2 powerful improvements exist for the measurement of the topological charge dis- 
tribution [TS] based on a generalization of the usual cluster algorithms 23 ,23]. The topological 
charge operator from [25] is used assigning an integer charge to each configuration even at finite 
lattice spacing. 

The continuum extrapolation of the cumulant differences will be done through 12 lattice 
spacings using the parameter values from |22) listed in table [T] The measured correlation 
lengths and topological susceptibilities are in agreement with those in [55] . In the present work 
O(10 8 ) configurations were generated at each volume and every 10 th was measured for the 
topological charge distribution and correlation length. The large number of configurations was 
necessary because there are huge cancellations between the various terms in the difference of 
cumulants, especially for Axa- The third difference, Axe, was already impossible to obtain 
with the current statistics. 

The results for the cumulant differences Ax2 and A\4 are shown on figure [T] Obtaining 
continuum estimates is not entirely trivial since the precise form of the leading and sub leading 
cut-off effects is not known a priori. Using the results of [H^HT] one may expect the leading 
corrections to be 0((a/L) 2 ) with possibly large logarithmic corrections. Fits of the form 

C+(a/Lf \^A 3 \og 3 {L/a)\ (12) 

with (n, m) — (0, 3), (1, 3), (2, 3), (0, 2) all work quite well with x 2 /dof values close to unity for 
A%2 and slightly higher, around 1.8 for Axa- The continuum extrapolated values agree in both 
cases among the four fit function choices and the four curves lie almost entirely on top of each 
other. In both cases the (n,m) = (0,2) choice is shown on the plots leading to continuum 
estimates C — 0.523(2) and 1.48(2) for L 2 A\2 and L 2 Ax4, respectively. Clearly, both values 
are non-zero. 
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L/a 


S/tt 


m(L)L 


L' 2 X2 


L' 2 X4 


L 2 X6 


L 2 A X 2 


L' 2 A X a 


60 


0.48490 


4.0017(14) 


1.2957(2) 


0.8812(8) 


-0.019(5) 


0.4145(8) 


1.069(5) 


80 


0.47260 


4.0032(19) 


1.4651(2) 


1.0292(8) 


-0.011(6) 


0.4359(7) 


1.143(6) 


100 


0.46370 


4.0007(19) 


1.6018(3) 


1.1512(9) 


-0.035(8) 


0.4507(9) 


1.186(7) 


120 


0.45680 


3.9939(20) 


1.7155(3) 


1.257(1) 


0.033(9) 


0.459(1) 


1.224(8) 


160 


0.44680 


4.0011(14) 


1.9214(4) 


1.444(1) 


0.16(1) 


0.477(1) 


1.28(1) 


200 


0.43950 


4.0015(17) 


2.0836(3) 


1.596(1) 


0.24(1) 


0.488(1) 


1.35(1) 


240 


0.43385 


3.9998(14) 


2.2208(3) 


1.729(1) 


0.40(1) 


0.492(1) 


1.33(1) 


320 


0.42545 


4.0010(17) 


2.4476(4) 


1.946(1) 


0.57(1) 


0.502(1) 


1.37(1) 


400 


0.41930 


3.9983(14) 


2.6259(4) 


2.118(2) 


0.71(2) 


0.508(2) 


1.41(2) 


480 


0.41455 


4.0014(19) 


2.7845(4) 


2.274(2) 


0.88(2) 


0.511(2) 


1.39(2) 


640 


0.40740 


4.0021(18) 


3.0347(4) 


2.521(2) 


1.07(3) 


0.514(2) 


1.45(3) 


800 


0.40210 


3.9952(19) 


3.2221(3) 


2.704(2) 


1.20(3) 


0.518(2) 


1.50(3) 



Table 1: Results for the first few cumulants and their differences for fixed physical volume 
m(L)L = 4. The bare parameters 6 are taken from [22] . 
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Figure 1: Continuum extrapolation for the first two cumulant differences multiplied by the 
volume, L 2 A\2 and L 2 A\4- 
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4 Summary and conclusion 



It has been known for a long time that the topological susceptibility in the two dimensional 
0(3) model is ill-defined in the continuum. Consequently the topological charge distribution 
P(Q) does not have a finite continuum limit. In this work it was determined precisely what 
part of P(Q) is actually divergent and what part of it is finite. The only divergent quantity is 
the first Fourier coefficient of the free energy density, 

h = ~ j /Wcos(tf)^, (13) 



while the remaining part JD n>1 (l — cos(w?))/ n is finite and non-zero. Hence the quantity 

IrW = fW - (1 - cob(0))/i (14) 

is finite and universal and one may consider the subtraction an additive renormalization. Sim- 
ilarly the renormalized partition function Zr({}) = exp(— V Jr^)) is finite and universal and 
related to the bare partition function by a multiplicative renormalization. Instead of subtracting 
/i it is sufficient to subtract only its divergent piece. The logarithmic singularity is expected 
to be volume independent Let us then denote this singular quantity by /i s . A natural renor- 
malization procedure is then the following: one defines the theory for non-zero $ by the action 

S{d) = S(d = 0) - idQ - (1 - cos{d))V f ls (15) 

and all resulting correlation functions related to topology (i.e. derivatives with respect to $) 
become finite. The last term in the full action above is a non-perturbatively generated counter 
term. It is important to note that the above renormalization does not mean that D itself gets 
renormalized, the bare d is still a physical quantity which does not require renormalization. It 
would of course be very interesting to check the volume independence of f\ s in lattice simula- 
tions. 

The finite quantities / n >i and Ax2m are not volume independent and are non-trivial func- 
tions of z = m(L)L. Since the model is exactly solvable at i? = it would be interesting to 
derive the first few cumulant differences Ax2m(z) from the exact solution or at least their value 
in the infinite volume limit. 

In any case the finite and non-zero cumulant differences naturally lead to the conclusion that 
# is a relevant coupling of the theory and the total topological charge operator Q is a relevant 
operator despite the ill-defined nature of the moments (Q 2m ). 

The original motivation was the study of a toy model mimicking confining, walking and 
conformal behavior in four dimensional gauge theories in order to study the severe systematic 
effects of the latter. It was proposed that increasing i? is analogous to increasing the number 
of flavors Nf because as $ goes from zero to ir the model goes from confining to walking and 
to conformal. In the toy model a suitable renormalized coupling is g\{L) = m(L)L which 
would then run with the finite volume L. A necessary condition for this analogy to hold was 
establishing precisely the divergence structure of the partition function at non-zero 

1 I thank Ferenc Niedermayer for pointing this out. 
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A particular difficulty of the gauge theory calculation can also be studied in the toy model. 
It is very difficult to distinguish numerically the following two cases: the theory with zero quark 
mass just below the conformal window and the theory with a small but non-zero quark mass 
just inside the conformal window. Both theories walk, the former for the usual reason of being 
just below the conformal window while the latter because even though it would be conformal 
for zero quark mass, the non-zero mass drives the coupling away from the would-be fixed point 
as soon as the running scale goes below the massive fermionic states. This phenomenon can be 
mimicked in the toy model by considering it at zero external magnetic field and $ = n — e and 
also at a small but non-zero external magnetic field and d = n. Both theories are expected to 
walk and it would be interesting to explore in the toy model what intrinsic features are different 
despite the similar behavior of the walking coupling constant. 

There are a couple of differences between the toy model and gauge theory though. Less 
important is the fact that while d does not enter the perturbative /3-function, Nf does. More 
significant is the fact that due to the i? — > — $ symmetry and periodicity by 2ir the two values 
fl — e and <& + e lead to the same continuum theory and it does not have an infra red fixed point 
(for non-zero e the coupling walks). This means that the zero of the /3-function at § = it is 
eliminated by arbitrary perturbations of 1? meaning that this zero is a second order zero, unlike 
in the gauge theory where generically the zero is expected to be first order and is preserved 
by small perturbations. Hence the 1) = it model is really analogous to a gauge theory which 
is exactly at the lower edge of the conformal window. It would be interesting to find a simple 
toy model which possesses all essential features and in addition the infra red fixed point is a 
first order zero of the /3-function and disappears by joining with a non-trivial UV fixed point 
as expected in gauge theory [28] . 
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